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ABSTRACT 

Using the published WMAP 5-year data, we first show how sensitive the WMAP power spec- 
tra are to the form of the WMAP beam. It is well known that the beam profile derived from 
observations of Jupiter is non-Gaussian and indeed extends, in the W band for example, well 
beyond its 12/6 FWHM core out to more than 1 degree in radius. This means that even though 
the core width corresponds to wavenumber I w 1800, the form of the beam still significantly 
affects the WMAP results even at I ps 200 which is the scale of the first acoustic peak. The 
difference between the beam convolved C/ and the final C/ is pa 70% at the scale of the first 
peak, rising to s» 400% at the scale of the second. 

New estimates of the Q, V and W-band beam profiles are then presented, based on a stack- 
ing analysis of the WMAP5 radio source catalogue and temperature maps. The radio sources 
show a significantly (3 — 4er) broader beam profile on scales of 10' — 30' than that found by the 
WMAP team whose beam analysis is based on measurements of Jupiter. Beyond these scales 
the beam profiles from the radio sources are too noisy to give useful information. Further- 
more, we find tentative evidence for a non-linear relation between WMAP and ATCA/IRAM 
95 GHz source fluxes. We discuss whether the wide beam profiles could be caused either by 
radio source extension or clustering and find that neither explanation is likely. We also argue 
against the possibility that Eddington bias is affecting our results. The reasons for the differ- 
ence between the radio source and the Jupiter beam profiles are therefore still unclear. If the 
radio source profiles were then used to define the WMAP beam, there could be a significant 
change in the amplitude and position of even the first acoustic peak. It is therefore important 
to identify the reasons for the differences between these two beam profile estimates. 

Key words: cosmic microwave background - cosmology: observations-early universe-space 
vehicles: instruments 



1 INTRODUCTION 

The WMAP satellite has produced some of the best support for the 
standard ACDM cosmological model. By measuring the first two 
acoustic peaks it has shown that the Univer se is spatially flat wit h 
tt A = 0.74 and H = 72 kms^Mpc -1 dfflnshaw et d]|2009h . 
The precision of the fit is impressive and rules out many co mpeting 
simple models suc h as the low Ho, fibaryon = 1 model of lShanka 
dl985ll2005tl2007t) . 

Of course, statistically precise measurements can also con- 
tain systematic errors which have to be guarded against. Such 
systematics include Galactic foregrounds which at the least cause 
mode coupling due to the incomplete sky (e.g. iHinshaw et al.l 
l2003l : IChonetalJl2004l) . There are also potentially more subtle 
syst ematics that arise from cosmological foregroun ds. For exam- 
ple, [MyersetaL I d2004h and lBielbv & Shanks) d2007t) have detected 
the SZ effect in the WM AP data by cro ss-correlating the CMB 
with rich cluster positions. Shanks! d2007» has also discussed the 
effect of foreground lensing, prompted by QSO lens ing results 
dMvers et alJl200ll2005l : lMountrichas & Shanksll2007l) . But SZ is 
unlikely to make a strong contribution to the first acoustic peak 
dHuffenberger et alj|2004) . Also lensing requires a high anti-bias 
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between galaxies and the mass distribution to have a significant ef- 
fect at the first peak which needs to be reco nciled with measures 
of bias from galaxy clustering dynamics (e.g. Ratclif fe"et al.ll 19981 : 
lHawkins et aU2003h . 

However, there are also many potential systematics involved 
with the WMAP instrument, although the WMAP team have taken 
care that the effects of such systematics are minimised. One ma- 
jor potential systematic concerns the question of the WMAP radio 
telescope beam profile. We shall see that even at the wavenumber 
/ « 220 of the first acoustic peak, the CMB power spectrum has 
significant dependence on the beam profile even in the highest res- 
olution W band. Here the W-band resolution quoted by the WMAP 
team is 12. '6 FWHM which is roughly equivalent to I « 1800. It 
is also noted that the beam is not Gaussian. Now the WMAP team 
have extensiv e papers devoted to the important question of measur - 
ing the beam dPage et al.l2003l : |jarosik et alj2007l : lHill et alj2009h . 
The standard method is to use their observations of bright sources 
such as the planet Jupiter to measure the beam profile. 

Here, after describing the WMAP5 data in [j2] we re-derive in 
ij3]the raw CMB power spectrum from the WMAP maps to show 
directly the effect of the beam. Then in ||4]we use radio sources to 
make new estimates of the WMAP beam and discuss other tests of 
the beam profile. In Sj5]we then make fits to the radio source beam 
profiles and use these to de-beam the WMAP5 data and explore the 
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Table 1. Summary of the WMAP sources listed as point sources in the 
Greenbank and PMN 5GHz catalogues. 



Band 


>l.Uy 


< l.Uy 


Total (> 2a) 


Q 


182 


165 


347 


V 


164 


153 


317 


W 


97 


84 


181 



range of power spectra that results. Possible reasons for the beam 
profile discrepancy and our conclusions are then presented in ij6] 

2 DATA 

2.1 WMAP5 maps and point source catalogue 

Here we use the five-year WMAP datasets which are available from 
the LAMBDA CMB website. The maps from the individual detec- 
tors in 5 frequency bands, K, Ka, Q, V and W are supplied. The 
FWHM of the 94 GHz W beam is 12. '6 compared to 19.'8 at V 
(61 GHz), 29. '4 at Q (41GHz), 37/2 at Ka (33GHz) and 49. '2 at 
K (23GHz). There are 10 differencing assemblies (DAs), namely 
Kl, Kal, Ql, Q2, VI, V2, Wl, W2, W3 and W4. The different DA 
maps can be cross-correlate d to obtain power sp ectra free of uncor- 
rected detector noise bias JHinshaw et aU2003l) . The Jupiter beam 
profiles for each DA and the correspon ding transfer funct ions are 
also given. The maps are in HEALPix (Gorsk i et alj|2005h format 
with iV s idc = 512 and 7V s idc = 1024. These give equal area pixels 
of dimension ~ 7' and ~ 3', respectively. 

We use the radio sources draw n from the WMAP5 point 
source catalogue dWrightetalJl2009t) . These sources have to be 
detected to > 5a in at least one WMAP band and their flux 
density is reported if they are detected at > 2<r in any of the 
other four WMAP bands. This gives a list of 390 sources to a 
limit of ~0. 5 Jy in each band . The source positions are accu- 
rate to ~ 4' dWright et al1 l2009). 365 out of 390 sources are pre- 
detected at 4.85 GHz in the Greenbank (GB6) northern sky survey 
|Gregorv et~ai1l 19961) an d the Parkes-MIT-NRAO (PMN) surveys 
(Griffith & Wright 1993). Here, we only use WMAP5 sources with 
4.85 GHz counterparts and exclude sources (12 out of 365) that 
were found to be resolved at ~ 4. '6 FWHM resolution of GB6 and 
PMN. Table 1 shows the number of these sources in each band, also 
split into those brighter or fainter than 1.1 Jy. 

From the optical identifications o f lTrushkirj J2003l) of the 208 
WMAP 1st year sources the survey contains 77 percent QSOs or 
BL Lac with the remainder being radio galaxies/ AGN. This is as 
expected given the dominance of flat-spectrum compact sources at 
the high WMAP frequencies. 

2.2 Ground-based 90-95GHz Radio Sources 

We shall compare WMAP W-ba nd fluxes with grou nd-based ra- 
dio source fluxes from ATCA dSadler et all 120081) and IRAM 
dSteppe et alii 19881) . The ATCA survey was made at 95 GHZ and 
the IRAM survey at 90GHz. The ATCA survey was based on 
sources selected at 20GHz. Of the 130 sources observed, 17 were 
detected at more than 2a by WMAP5. The IRAM survey ob- 
served 294 sources at 90 and 230 GHz, targetting sources which 
are brighter than 1 Jy at 5 GHz. Here 66 sources were detected at 
more than 2a by WMAP5. At these high frequencies the sources 
are mainly QSOs, BL Lacs or blazars. Many of the sources in the 
ATCA and IRAM surveys are variable and so where this is an issue 
we shall use the average source fluxes in our comparison with the 
WMAP fluxes. 
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Figure 1. The beam convolved C; measured from the Wl and W2 maps 
(green line) compared to the standard de-beamed WMAP result as pre- 
sented by Nolta et al. (blue line). Also shown are C; measurements of a 
full-sky CMB simulation (WMAP5 best-fit ACDM) smoothed with a 12/6 
FWHM Gaussian beam (black line), and a similar simulation now smoothed 
with the Wl and W2 Jupiter beams (red line). Although the beam convolved 
C; (green) and the Jupiter beam-smoothed standard simulation C; (red) 
agree, the difference between these and the Nolta et al result (blue) shows 
the large effect of debeaming even at the scale of the first acoustic peak. 

3 DERIVING THE BEAM CONVOLVED Ci 

We now analyse the WMAP W-band data to make an initial esti- 
mate of the power spectrum from the W band. To reduce the ef- 
fect of correlated detector noise which would result from an auto- 
correlation of an individual CMB map, we make a cross-correlation 
of the maps from indep endent detectors. W e derive the result by 
using the PolSpice code (jSzapudi et al j200ll) to cross-correlate the 
Wl and W2 maps. The default WMA P5 temperature power spec- 
trum mask (KQ85, iNolta et al.| [2009) is used and the cut-sky cor- 
rected angular pow er spectrum, C[, is obtained directly from the 
PolSpice code (see IChon et alj|2004[) . Hereafter, we shall call an 
angular power spectrum after a cut-sky and pixel transfer function 
(see Eq. |2j correction a 'beam convolved Ci . In Fig. [T] we im- 
mediately see that the beam convolved C; (green line) is not only 
drastically smoothed at the position of the second and third peaks 
but there is also a significant effect at the position of the first peak 
at I « 220 in that the amplitude of the standard ACDM result (blue 
line) is ~ 70% higher. The reason for this is seen in Fig. |2f where 
the beam profile from the Jupiter observations using the Wl detec- 
tor are shown. It can be seen that the beam is not Gaussian and has 
a 9~ 3 power-law tail o ut to > 1°. 

|PageetaI] ( l2003h give the relation between the beam transfer 
function, bi, and the radial beam profile, b s (6), as, 

bi = 2n J b s {6)Pi{cos9)dcos9/£l B (1) 

The de-beamed cross-power spectrum measured from DA i and j 
is then given from the measured G{ by 

d = C'Jblbfpl (2) 

where pi is the pixel transfer function supplied with the HEALPix 
package. For iV s id c = 512, the pixelisation lowers the measured 
C'i by ~ 1 and 10 percent at I 200 and 500, respectively. 
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Figure 2. (a), (b) and (c); The raw radio source beam profiles for the WMAP sources compared to the beam profiles from Jupiter (solid lines) and random 
realisations (dotted lines) for Ql, VI and Wl. (d), (e) and (f); The zero-offset subtracted (see text) and normalised beam profiles for Ql, VI and Wl. The 
radio source profiles for all compact WMAP sources are shown as diamonds. Profiles derived from WMAP sources with flux brighter (fainter) than 1.1 Jy are 
shown as asterisks (stars). A Gaussian is shown as a dot-dashed line and empirical fits to the radio source profiles are shown as orange and green dashed lines. 



If we use Eq.[2]with the Jupiter beam transfer function from 
the WMAP team, we find that we get back to the usual ACDM 
model (green and orange lines in Fig. [4]>. The black line shows the 
Ci measured from a full-sky CMB simulation (WMAP5 best-fit 
ACD M model) after sm oothing by a Gaussian beam, using synfast 
code dGorski et al.ll2005t) . The red line shows the effect of the sim- 
ilar simulation now smoothed with W1W2 Jupiter beam profiles. 
The latter shows excellent agreement with our beam convolved Ci 
measured from the Wl and W2 maps. The W1W2 Ci is noisier 
than the simulation due to radiometer noise in the data. The effect 
of the Jupiter beam compared to the Gaussian beam is thus very 
significant in decreasing the height of the first peak. We also see 
that when the Jupiter beam is used, the ACDM model does give an 
accurate fit to the beam-convolved Ci . Thus when we use the same 
parameters as the WMAP team, we reproduce the WMAP result. 



4 TESTING THE WMAP BEAM PROFILE 
4.1 Beam profiles via point sources 

We then estimated a beam from the radio sources by making 
a stacking analysis of WMAP5 temperature maps around radio 
source positions. The extended foreground emission regions are 
excluded fr om the temperatur e maps using the 'point source cata- 
log mask' JWright et al j2009t) . We calculated the average AT (per 
49 arcmin 2 pixel) in annuli as a function of angular distance, 8, 
between radio source position and the pixel centre. In the first in- 
stance we show the raw cross-correlation function for the Q, V and 
W bands in Fig. [2] split into bright (^ 1.1 Jy) and faint (< 1.1 Jy) 
WMAP5 source sub-samples. The errors on the radio source pro- 
files are jackknife errors. These are estimated using six equal-area 
sub-fields, given by a 2 ]K (6) = (Nj - 1)([AT(0) - AT(6>)] 2 ), 



where Nj — 6 and AT is the stacked temperature measured from 
five out of the six sub-fields. 

We see that the fainter source profiles appear to agree with 
the brighter source profiles at scales of 9 m 30' but have signif- 
icantly lower peak values. This is most clearly shown in the un- 
renormalised profiles shown for the bright and faint Q,V, W band 
sources in Fig. [2^, b and c. Although noise may be an issue for 
the faintest sources, this suggests that there may possibly be some 
form of non-linearity in the WMAP beam. We also note that the 
profiles from both bright and faint sources show a positive offset at 
the 0.01-0.02 mK level. The offset shows an increasing trend from 
1° — 5°. The main uncertainty in estimating WMAP beam profiles 
from these data is in subtracting this offset at scales > 1°. 

Since WMAP has significant sidelobes stretching to m 90° 
( iBarnes et alj2003h , there was a possibility that the offsets are also 
part of the beam. However, when we distributed points at random 
in the masked region and used these as our centres for our stacking 
analysis we also found a similar offset (dotted lines in Fig. [2^, b 
and c). This makes it look like the offset is not associated with 
the existence of sources and hence not associated with the WMAP 
beam. Our Monte Carlo simulations (see below) shows that these 
offsets are caused by the CMB fluctuations. We therefore employed 
a 'photometry' approach for the stacking analysis where we have 
subtracted the WMAP flux in an annulus at 1° < < 2° for the W 
band and proportionately bigger annuli in the V and Q bands. Using 
sky annuli close to the sources will clearly improve background 
subtraction in the presence of local background fluctuations. 

The resulting WMAP radio source beam profiles for Q, V and 
W bands are shown in Fig. |2jl, e and f. The profiles have been 
renormalised (~ 10% statistical uncertainties in the normalising 
factors) to fit the peak in the WMAP Jupiter beam profile at 6 < 4' 
and this profile is also shown. For each band we also compare the 
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profiles to a Gaussian beam with the FWHM as indicated in the 
plot. We see that on average the profiles from the radio sources 
are broader than the Jupiter profile in the W, V and Q bands. In 
the lower frequency, lower resolution K and Ka bands the radio 
source profiles fit the Jupiter beam better, indeed almost perfectly 
(not shown here). Clearly, given the size of the errors there is little 
information from the radio sources on the beam profile at 9 > 30'. 
Fig.|2]again shows the WMAP radio sources divided into faint and 
bright sources, split at 1.1 Jy. In the W and V bands particularly 
we again note that the fainter sources appear to be wider than the 
brighter sources. We also find similar results for W2, W3, W4, V2 
and Q2 but choose not to include them here for clarity. These de- 
viations from the WMAP Jupiter beam are puzzling and we now 
check to see if they could be caused by systematics. 

At the refe ree's request , we made 100 Monte Carlo simula- 
tions following Iwi-ight et alj J2009h to check our results. In sum- 
mary, simulated maps are constructed to include point sources sam- 
pled from a power law N(> S) distribution with spectral charac- 
teristics as seen in the data. The temperature map for each band is 
then smoothed with the WMAP Jupiter beam profile before being 
added to a simulated CMB map including radiometer noise. We 
then applied the five-band detect ion method following procedures 
described bv lWright et alj {2009). Applying our stacking analysis 
described above, we found that even profiles as narrow as the W- 
band Jupiter profile can be accurately retrieved. The flux depen- 
dence of measured profiles were small with only a hint of possible 
Eddington bias in the faintest bin. The pixelisation effect is also too 
small to explain the wider profiles seen here. Furthermore, the esti- 
mated uncertainties using these simulations are consist ent with the 
Jackkn ife error estimates. Further details are given bv lSawangwij 
J2010h who, further, finds similarly wide profiles for NVSS radio 
sources flux limited at 1.4 GHz, where any Eddington bias would 
be negligible. 

We have checked the likely contribution of radio source clus- 
tering to the beam pro files, using the clus tering analysis of the 
NVSS radio survey by lOverzier et alj d2003h . At S > 200 mjy 
where the sky density of NVSS sources is n ~ 0.6 deg -2 , 

w{9) = 3 x 10~ 5 6r 3 ' 4 + 6.6 x 10" 3 6» -0 ' 8 . (3) 

This is a 2-power-law form which changes slope at 9 ~ 0.1°. 
At smaller scales, double-lobed radio sources split into two com- 
ponents dominate while at larger scales source-source clustering 
dominates. We first calculate the excess number of sources in an 
annulus of area AA at radius 9 from an average source, N ex (9) = 
w(9)nAA. The excess flux/temperature per unit area in the pro- 
file in the annulus is then given by AT ex = N ex x //A A where 
/ is the average source flux. For a Gaussian point source of cen- 
tral intensity/temperature per unit area, To, and width, a, the flux 
is 2na 2 T . Therefore in this case, AT ex (9) = w(9)2nna 2 T . We 
find AT ex « 3 x 10~ 4 Tb which is a negligible contribution in 
explaining the excess in Fig. [2f at this scale, if Tb is taken to be 
the central profile va lue. Taking the parameters for 100 mjy from 
lOverzier et al. (2003) makes the effect even smaller. These results 
are also likely to be upper limits for the WMAP sources which only 
have a density of n ~ 0.01 deg~ 2 and an average 95 GHz flux of 
500 mJy. We conclude that radio source clustering is not likely to 
be an issue for our radio source beam profiles. 

We conclude that in the W and V bands and probably the Q 
band, the average radio source profile is wider than the Jupiter beam 
and the fainter sources may show a wider profile than the brighter 
sources. For Wl and S > 1.1 Jy sources, the beam profile mea- 
sured here rejects that of Jupiter with 4.0, 3.0 and 3.5ct significance 
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Figure 3. The comparison of the ATCA (diamonds) and IRAM (asterisks) 
source flux densities with the WMAP W-band data. 

for 9 = 12. '6, 20' and 31/6. These become 4.4, 3.2 and 2.8a when 
Monte Carlo errors are used instead. Note that the pixelisation has 
been taken into account when estimating these significances. 

4.2 Comparison with ground-based fluxes 

W e now make a check of the WMAP5 W band fluxes as presented 
bv lWright et alj {2009) in their Table 1. We checked these against 
the ATCA and IRAM source flux densities. The comparison in Fig. 
[3]shows that for both surveys, the brighter sources with fluxes > 3 
Jy are about a factor of 1.5 fainter in the WMAP source list than in 
the ATCA or IRAM lists. The agreement between the ATCA and 
IRAM fluxes appears better than for WMAP, if we use WMAP as 
an intermediary between these two surveys. If the scale error is due 
to WMAP, then this might suggest that there is a non-linearity in 
the WMAP flux scale. It could mean that a narrower WMAP beam 
at brighter fluxes is missing a significant amount of flux in the tail 
of the beam profile. 

5 IMPACT ON THE DE-BEAMED Ci 

Finally, we use the information from our radio source beam profiles 
to judge what the effect might be on the de-beamed WMAP Ci. 
Unfortunately we will have to extrapolate our radio source fits in 
the regime beyond ~ 1° out to 5° because of the large errors on the 
radio source beam profile in this range. These results can therefore 
only be used to indicate the sensitivity of the C; to the beam profile 
and should not be regarded as alternative Ci estimates. We first 
make an extrapolation where we fit the small-scale (total sample) 
beam profile points and then extrapolate continuing with the power- 
law as shown by the green dashed lines in Fig. [2p, e and f. We 
also made a more conservative extrapolation where we again fit 
the small-scale data but then extrapolate continuing parallel to the 
Jupiter beam profiles at large scales (orange dashed lines). 

The range of the radio source de-beamed Ci is shown by the 
two red lines and two cyan lines (for W and V bands) in Fig. [4] The 
most conservative profile model is ~ 50% higher than the Jupiter 
de-beamed Ci (green and orange lines) at the scale of the first peak. 
But the most extreme model is now a factor of 2-3 higher even at 
I — 220 than the standard model power spectrum. We note that 
it has been possible to derive consistent Ci's between the V, W 
(and Q) bands, although we accept that this is due to the freedom 
we have in extrapolating our radio source beam profiles beyond 
9 ~ 30'. It seems that if the radio sources are indicating a wider 
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Figure 4. The de-beamed CMB power spectra. The blue diamonds show 
the WMAP team's result. The cross-Q for the Wl and W2 (VI and V2) 
maps and de-beaming with the Jupiter beam profile is shown as the green 
(orange) line. The same result but now de-beamed using the 12. '6 FWHM 
Gaussian is shown as the magenta line, significantly different from the result 
using the Jupiter beam profile to de-beam. The same results but now using 
the profiles from the compact radio sources are shown as the red and cyan 
lines for W and V band respectively. The differences in the two sets of de- 
beamed W and V spectra are due to the difference in extrapolations of the 
radio source beam profiles beyond 9 = 30' (see Fig. [2^ and[2j.) 

beam profile, then the systematic uncertainty in the beam at the 
largest scales will dominate the error budget of the Ci even at the 
scale of the first acoustic peak. These larger errors would then allow 
a wider range of cosmological models to be fi tted, including models 
where the first peak lies at / as high as 330 dShanks & Saw angwil 
HH3). 

6 DISCUSSION AND CONCLUSIONS 

Clearly it is important to understand why the radio source profiles 
are so wide in the Q, V and W bands. If there is a correlation be- 
tween beam width and source flux then it will be wrong to use 
Jupiter to debeam the CMB power spectra because in the W band, 
for example, the ~1 Jy radio sources are much closer to the ~ 0.5 
Jy rms flux of the CMB fluctuations than the 1200 Jy flux of Jupiter. 

The non-linearity shown by the WMAP source fluxes com- 
pared to independently measured ATCA/IRAM fluxes is supporting 
evidence of non-linearity in the WMAP data calibration. It is pos- 
sible that somehow the variability of the radio sources at W have 
caused problems that would not apply to the CMB. In passing, we 
note that the smaller than ex pected SZ decrements from WMAP 
obser vations of rich clusters dMvers et alj|2004l ; iBielbv & Shanks! 
1 20071) may also be explained by a wider than expected WMAP 
beam at W. If so, this would argue that the non-linearity affects 
variable and non-variable WMAP sources alike. 

In considering possible causes of WMAP non-linearity, we 
first note that detector saturation is unlikely to be the problem since 
this would lead to the brighter sources having a wider profile than 
the fainter sources, which is opposite to what is observed. However, 
Jupiter, being a moving source, has to be dealt with in a different 
way to the radio sources and the CMB fluctuations in the maps. 
This means that if there was a problem in the WMAP analysis, it 
would be necessary to check any filtering that is done to the maps. 
Otherwise, we do not understand the reason for the difference be- 
tween the Jupiter and radio source beam profiles. 



We conclude that; 

• The WMAP power spectrum is heavily dependent on the beam 
profile. Indeed even the first acoustic peak at I S3 220 is very de- 
pendent on the form of the profile at 1° — 2° where the profile is 
only S3 0.1 percent of its peak value. 

• The radio point sources detected by WMAP in the Q, V and W 
bands generally show a broader beam profile than the Jupiter beam 
used by the WMAP team. For example, using bright point sources, 
our Wl beam profile rejects the Jupiter beam with > 99.5% confi- 
dence. 

• There may be evidence for a flux dependent effect within the 
WMAP data in that fainter radio sources appear to have systemati- 
cally broader profiles than brighter sources, although the faint data 
are noisy. 

• Non-linearity in the WMAP flux scale may also be indicated 
by comparisons of WMAP radio source fluxes with ATCA and 
IRAM fluxes which show 50 per cent reduced flux from WMAP. 

• Further arguments against possible systematics such as Ed- 
dington bias affecting our results come from simulation checks and 
NVSS source samples selected at frequency where CMB fluctua- 
tions are subdominant (see Sawangwit 2010). 

• The systematic errors on the WMAP C; due to the beam may 
be much larger than previously expected and in turn, this means 
that the systematic error on the best fit cosmological model may 
also be larger. It will be interesting to see if a revised estimate of 
the WMAP beam profile then allows a simpler cosmological model 
to be fitted than ACDM. 
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